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ABSTRACT 


In this study we have develop a basic deterministic mathematical model to 
investigate SEIR Model and Simulation for controlling malaria Diseases 
Transmission without Intervention Strategies. The model has seven non¬ 
linear differential equations which describe the spread of malaria with 
three state variables for mosquitoes populations and four state variables 
for humans population and to introduce the model without intervention 
strategies. 

The models are analyzed qualitatively to determine criteria for control of a 
malaria transmission, and are used to calculate the basic reproduction R_0. 
The equilibria of malaria models are determined. In addition to having a 
disease-free equilibrium, which is locally asymptotically stable when the 
R_0<1, the basic malaria model manifest one's possession of (a quality of) 
the phenomenon of backward bifurcation where a stable disease-free 
equilibrium co-exists(at the same time) with a stable endemic equilibrium 
for a certain range of associated reproduction number less than one. The 
results also designing the effects of some model parameters, the infection 
rate and biting rate. The numerical analysis and numerical simulation 
results of the model suggested that the most effective strategies for 
controlling or eradicating the spread of malaria were suggest to use 
insecticide treated bed nets, indoor residual spraying, prompt effective 
diagnosis and treatment of infected individuals. 

KEYWORDS: Malaria, Basic reproduction number, Stability analysis, 
Existence of Backward bifurcation analysis, Endemic equilibrium point 


1. INTRODUCTION 

Malaria is an infectious disease and is life threatening for 
human beings worldwide. Parasite is an organism that 
lives on or inside a human body from which it gets its food. 
Malaria is caused due to a parasite called Plasmodium. 
Plasmodium parasite is transmitted into human body 
when an infected female anopheles mosquito makes bites. 
Plasmodium parasites making the human liver as their 
home multiply their population and start infecting red 
blood cells of the human. A variety of plasmodium 
parasites exist. Mainly four types of plasmodium cause 
malaria disease among the human viz., falciparum, vivax, 
ovale and plasmodium malaria [60], 

Malaria is an infectious disease and is life threatening for 
human beings having a huge social, economic, and health 
burden. Malaria transmission occurs in all over the 
worldwide. Globally, an estimated 3.2 billion people are at 
risk of being infected with malaria and developing disease, 
and 1.2 billion are at high risk (>1 in 1000 chance of 
getting malaria in a year). According to the latest 
estimates, 198 million cases of malaria occurred globally 
in 2013 (uncertainty range 124-283 million) and the 
disease led to 584,000 deaths (uncertainty range 367,000- 
755,000).The burden is heaviest in the WHO African 
Region, where an estimated 90% of all malaria deaths 
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occur, and in children aged under 5 years, WHO account 
for 78% of all deaths [45], 

In our country Ethiopia is a major public health problem 
and has been reported the last five years (2002-2008) the 
proportion of malaria in outpatient department, admission 
and in-patient deaths has been increasing with the highest 
being recorded in 2003 and 2004. In 2008 malaria was 
still the first leading cause of health problem accounting 
for 48% of outpatient consultations, 20% admissions and 
24.9% inpatient deaths. According to FMOH reports, 
approximately 70,000 people die of malaria each year in 
Ethiopia [10], 

Malaria is a life threatening infectious disease caused by a 
parasite called Plasmodium which is transmitted through 
the bites of infected female anopheles mosquitoes. There 
are four different species causing the human malaria 
disease plasmodium falciparum, plasmodium vivax, 
plasmodium ovale and plasmodium malaria ([44],[50]). 
The plasmodium parasite is injected into the human 
bloodstream in the form or stage or life cycle known as 
sporozoite. The parasites go through a complex life cycle 
inside the hosting human body and they live at various 
stages both in liver and red blood cells. From time to time 
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the parasites pass through various stages of their life cycle 
and during which numerous human red blood cells are 
destroyed. From the listed four plasmodium parasites in 
our country Ethiopia are plasmodium falciparum and 
plasmodium vivax, accounting for 60% and 40% cases 
respectively present in country Ethiopia widely [11], 

At this stage the disease generates and develops its 
symptoms in the infected human body. Eventually, the 
parasites become gametocytes which are in turn taken by 
mosquitoes that bite the human host. Inside the mosquito, 
the gametocytes mature, reproduce sexually, and migrate 
into the mosquito's salivary glands, at which stage the life 
cycle is repeated. For some species of Plasmodium, the 
parasites may persist in the liver for months or years, 
resulting in chronic and recurring eruptions of merozoites 
that correspond to episodes of fever and sickness. 

The effect of malaria disease varies with the infecting 
variety species of Plasmodium and also with prior health 
and immune status of the individual. Typically malaria 
disease causes fever and chills together with headaches, 
vomiting and diarrhea. It may also cause long-term 
anemia, liver damage and neurological damage. The most 
dangerous falciparum parasite can cause cerebral malaria 
which causes frequently a fatal condition involving 
damaging the brain and central nervous system. The 
survived people from the cerebral malaria may too 
experience brain damage. 

Now a days, although malaria deaths do not occur as often 
as previously, but still it remains a major public health 
problem and it is too early to reach any firm conclusion 
about the possibility of achieving MDGs, because of 
resistance of the parasite to antimalarial drugs, the 
complexity of disease, expensiveness of the control 
program, seasonal variability nature of the disease [44], 

In the recent time, significant resources and control 
programs have been made available worldwide. The aim is 
to reduce malaria infected cases and prevalence or gain 
upper hand over the disease. Different strategies and 
programs with varying effectiveness and efficiency are 
being adopted to control malaria disease. Comparative 
knowledge of these existing programs is necessary to 
design and organize any new and useful and cost effective 
procedure to control malaria epidemic ([44],[52]]. 

The National Strategic Plan for Malaria Control and 
Prevention in Ethiopia (NSP] 2006-2010 aimed to rapidly 
scale-up malaria control interventions to achieve a 50% 
reduction of the malaria burden, in line with global Roll 
Back Malaria (RBM] [35] partnership objectives. The 
status of coverage of the major interventions was 
measured in the Malaria Indicator Survey [MIS] 2007. The 
MIS 2007 results show tremendous achievements by 
Ethiopia’s malaria control program. Thus, between 2005 
and 2007, insecticide-treated net (ITN] coverage 
increased 15 fold, with ITN use by children under five 
years of age and pregnant women increasing to nearly 
45% in malaria-endemic areas and to over 60% in 
households that owned at least one ITN. Overall, 68% of 
households in malaria-endemic areas were protected by at 
least one ITN and/or indoor residual spraying of 
households with insecticide [IRS], It is believed that the 


vector control interventions have contributed greatly to a 
reduction in the burden of the disease. More than 20 
million LLINs have been distributed to 10 million 
households between 2005 and 2007. With respect to IRS 
activities, evidence shows that 30% of IRS-targeted areas 
were sprayed in 2007 and in 2008 the coverage increased 
to 50%. So far, the main vector control activities 
implemented in Ethiopia include IRS, LLINs and mosquito 
larval source reduction. 

The Malaria Vector Control Guidelines also addresses 
vector control interventions found to be effective in past 
decades. The insecticides commonly used in the country 
include dichloro-diphenyl-trichloroethane (DDT], 
Malathion and deltamethrin. Due to resistance of malaria 
vectors to DDT, the use of this Insecticide for IRS has been 
discontinued in 2009. Deltamethrin is currently being 
used as an interim substitute insecticide for DDT in IRS 
operations. However, the selection of insecticides for IRS 
use in Ethiopia will be determined annually based on the 
insecticide resistance pattern of the vectors and other 
factors. Environmental management, supported by active 
participation of the community and use of larvicides are 
other preventive measures addressed in this guideline. 
The guideline incorporates the three major vector control 
measures, namely environmental management, IRS, and 
LLINs [11], 

Efforts to reduce malaria transmission have led to the 
development of efficient vector control interventions, 
particularly insecticide treated nets (ITNs], indoor 
residual spraying(IRS],and larvicide ([53],[54]]. The ITNs 
include conventional nets treated with a WHO 
recommended insecticide and long-lasting insecticidal 
nets. Note that larva is an immature form of an insect and 
larvicide is a chemical used to kill larvae. These 
interventions are used in malaria endemic countries 
especially those in sub-Saharan Africa and have led to 
reduction in malaria morbidity and mortality 
substantially. However, malaria epidemic continues to 
claim hundreds of thousands of lives every year, thus 
necessitating a continued control effort to fight against the 
disease ([55], [56]]. 

Malaria has been considered as a global issue. 
Epidemiologists together with other scientists invest their 
efforts to understand the dynamics of malaria and to 
control transmission of the disease. From interactions 
with these scientists, mathematicians have developed 
tools called mathematical models. These models are used 
significantly and effectively for giving an insight into the 
interaction between the humans and mosquito population, 
the dynamics of malaria disease, control mechanisms of 
malaria transmission and effectiveness of eradication 
techniques. 

Mathematical models are particularly helpful as they 
consider and include the relative effects of various 
sociological, biological and environmental factors on the 
spread of disease. The models have played a very 
important role in the development of malaria 
epidemiology. Analysis of mathematical models is 
important because they help in understanding the present 
and future spreads of malaria so that suitable control 
techniques can be adopted. 
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The SEIR is a simplest mathematical model and has four 
classes or compartments Susceptible, Exposed, Infected 
and Recovered. The effect of controlling technique in the 
spread of malaria is analyzed. Using these notations, eight 
classes of compartmental models are possible, SI, SIS, SEI, 
SEIS, SIR, SIRS, SEIR and SEIRS ([14],[29],[57]]. For 
example, in an SEIRS model, a fraction of the susceptible 
(S) population gets exposed (E) to infection, a part of 
which then becomes infectious (I], Some from the I class 
recover from the disease, and become part of the R class 
with temporary immunity. When immunity is lost, they 
become susceptible to pathogen attack again, and enter 
the S class. The simulation studies of the model with 
variable values of sensitive parameter of the spread of 
malaria are performed and the results are incorporated. 
The necessary conclusions have been drawn. 

2. The SEIR Malaria Model without Intervention 
Strategies 

2.1. The model formulation 

In this study we formulate that of similar model [58] and 
[60] describing the transmission of malaria. The malaria 
model divides the human population into four classes and 
with assumptions about the nature and time rate of 
transfer from one classes to another. We consider the total 
population sizes denoted by N h (t) and N m (t ) for the 
human hosts and female mosquitoes, respectively. We will 
use the SEIRS framework to describe a disease with 
temporary immunity on recovery from infection. SEIRS 
model indicates that the passage of individuals is from the 
susceptible class, S, to the exposed class, E, then to 
infective class ,7, and finally to the recovery class, R. S(t) 
contains humans those do not have malaria disease but 
are likely to be bitten by infected female anopheles 
mosquitoes causing malaria parasite at time t, or those 
susceptible to the disease. 

Many diseases like malaria have what is termed a latent or 
exposed phase, E(t), during which an individual is said to 
be infected but not infectious. 7(t) contains humans those 
are already infected and got malaria disease. People come 
into infected class from susceptible class. This is done 
through infecting the susceptible mosquitoes. The 
dynamic transmission of the malaria parasite between and 
among individuals in both species is driven by the 
mosquito biting habit of the humans. R(t) contains the 
people who recover from the malaria disease and return 
to normal status of health or individuals who have 
recovered from the disease. These humans cannot 
transmit the infection to mosquitoes as we assume that 
they have no plasmodium parasites in their bodies. 

The transfer rates between the subclasses are composed 
of several epidemiological parameters, we explained that 
susceptible human bitten by an infectious anopheles 
mosquito may become infected with a finite probability 
that depends on the abundance of infectious mosquitoes 
and human populations [20], The model assumes straight 
forward to the occurrence of the same kind that 
susceptible individuals get infected through biting with 


infected mosquitoes. The susceptible human population is 
increased by recruitment (birth and immigration] at a 
constant rate y. 

All the recruited individuals are assumed to be naive when 
they join the community. Infected immigrants are not 
included because we assume that most people who are 
sick will not travel. When an infectious female anopheles 
mosquito bites a susceptible human, there is some 
probability of transmission of infection from an infectious 
mosquito to a susceptible human d mh . 

The parasite then moves to the liver where it develops 
into its next life stage. The infected person will move to the 
exposed class. After a certain period of time, the parasite 
(in the form of merozoites] enters the blood stream, 
habitually the indication of malaria disease recognizable 
when the symptoms begin on the human body. Then the 
exposed individuals become infectious and progress to 
infected state at a constant rate fi m . We exclude the direct 
infectious-to-susceptible recovery by assuming that the 
individuals do not recover by natural immunity. This is a 
true life to simplifying assumption because most people 
have some period of immunity before becoming 
susceptible again. After some time, individuals who have 
experienced infection may recovered with natural 
immunity at a constant rate t and move to the recovered 
class. The recovered individuals have some immunity to 
the disease and do not get clinically ill. Since disease- 
induced immunity due to malaria is temporary a fraction 
c|> of individuals leave the recovered state to the 
susceptible state. We make the simplifying assumptions 
that there is no immigration of the recovered humans. 
Humans leave the population through natural death m and 
the infected humans have an additional disease-induced 
death rate constant 8. The disease-induced rate is very 
small in equivalent with the recovery rate. 

We divide the mosquito population into three classes: 
susceptible S m exposed E m and infectious I m . Female 
anopheles mosquitoes (male anopheles mosquito is not 
included in the model because only female mosquito bites 
humans for blood meals] enter the susceptible class 
through natural birth at a rate p. Susceptible mosquitoes 
become infected by biting infectious humans at a 0. The 
parasites (in the form of gametocytes] enter the mosquito 
with probability 0 hm , when the mosquito bites an 
infectious human, and the mosquito moves from the 
susceptible to the exposed class. After some period of time, 
dependent on the ambient temperature and humidity, the 
parasite develops into sporozoites and enters the 
mosquito's salivary glands, and the mosquito progresses 
at rate (3 m , from the exposed class to the infectious class. 
We assume that the infective period of the vector ends 
with its death, and therefore the vector does not recovered 
from being infective [1], The mosquitoes leave the 
population through natural death. Its caused by natural 
death rate and insecticides is or. The rate of infection of 
susceptible individual is a h , and the rate of infecting a 
susceptible mosquito is a m . 
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Table 1 Variables of the basic malaria model 


Variables 

Description 

s h (t) 

Number of humans insusceptible compartment at time t 

E h (t) 

Number of humans in exposed class at time t 

4 ft) 

Number of humans in infected compartment at time t 

R(t ) 

Number of humans in recovered compartment at time t 

S m (t) 

Number of mosquitoes in susceptible compartment at time t 

Emit) 

Number of mosquitoes in exposed class at time t 

4j (t) 

Number of mosquitoes in infected class at time t 

N h (t) 

Total human population at time t 

NmXt ) 

Total mosquito population at time t 


Table 2 Parameters and their interpretations for the malaria model 


Parameter 

Description 

V 

Natural birth rate of humans 

P 

Natural birth rate of mosquitoes 

a h 

Transfer rate of humans from susceptible to infected compartment 


Transfer rate of mosquitoes from susceptible to infected compartment 

p 

Natural death rate for humans 

s 

Death rate of humans due to disease-induced 

Cl) 

Death of mosquitoes caused by natural death rate and insecticides 

P, 

Transfer rate of humans from the exposed class to the infected class 

t 

Transfer rate of humans from Infected to recovered class 

(p 

Transfer rate of humans from recovered to susceptible compartment 

Qm.h 

Probability of transmission of infection from an infectious mosquito to a susceptible human 

8h.m. 

Probability of transmission of infection from an infectious human to a susceptible mosquito 

Pm 

Transfer rate of mosquitoes from the exposed class to the infected class 

0 

Susceptible mosquitoes bite infected humans with this rate. 

Also, infected mosquitoes bites susceptible humans with the same rate 


Figure 1 The flow chart for transmission of malaria disease 

HUMAN POPULATION MOSQUITO POPULATION 



2.2. Mathematical formulation of SEIR model 

Applying the assumptions, definitions of compartmental 
variables and parameters described in tables 1 and 2, the 
system of non-linear differential equations which describe 
the dynamics of malaria transmission with controlling 
measures are formulated and presented in this section. 

dS h /dt — y/+ <pR — iiS h — a h S h 

dEfr j dt — oCfrSfr -(- fifrEfr [tEfr 

dljdt = p h E h - Tl ft - (p + $ft)Ift 

dR/dt = Tl h — cpR — pR (1) 

dS m f dt — p — cc m S m — 

d-E m / dt — <x m S m + p m E m — coE m 

d\ m /dt p m S m mi m 


The initial conditions of the system of equations (1) are 
given by Sft(O) = S h0 ,E h ( 0) = E hQ , 4(0) = I h0 ,R(0) = R 0 
- ^m(O) = S m0 , E m ( 0) = E mQ and / m (0) = I m0 . Also, we 
have used in equation (1) that a h = (O mh 0 \ m /N h ) 
and a m = ( 9 hm 0 l h /N h ). The term a h denotes the rate at 
which the susceptible humans become infected by 
infectious female mosquitoes. Similarly, the term a m 
denotes the rate at which the susceptible mosquitoes 
become infected by infectious humans. The rate of 
infection propagated to susceptible humans by infected 
mosquitoes is dependent on the total number of humans. 
Similarly, the rate of infection propagated to susceptible 
mosquitoes by infected humans is dependent on the total 
number of humans ([18], [60]). 
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2.3. Analysis of SEIR model (Analysis of the Model 
without intervention strategies) 

We now analyze the SEIR model in order to show the two 
controlling methods considered here have substantial 
impact on controlling the transmission dynamics of 
malaria disease. In fact, the disease will be completely 
eradicated if the controlling methods are implemented 
effectively. The two controlling mechanisms proposed 
here have such a big potential. We consider now the 
solutions of the system of non-linear differential equation 
(1). We understand that the interpretations of these 
solutions must be biologically meaningful. Hence it is easy 
to identify that the feasible region of system (1) is M+. The 
seven dimensional solution space shows that all the 
solutions are positive. Hence, the feasible region 
containing all the solutions of the system of equations (1) 
is given by the set/2 = {( S h ,E h , I h , R,S m , E m , I m ) £ IR+}. 
Here the quantities S h ,E h , I h ,R,S m ,E m ,l m are all non¬ 
negatives. Further the total human and mosquito 
populations are represented by N h and Al )n they have the 
upper asymptotic values ( y//p ) and (p/a>) respectively. 
Therefore, the region /2 is positively invariant i.e. solutions 
remain positive for all the temporal values. Thus, the 
model (1) is biologically meaningful and mathematical 
well-posed or well present in the domain 0. 

On summing up all the individual equations from (1) of 
the system (1), it is straight forward to get ( dN h /dt ) = 
(i//- pN h - 8 h I h ). Here the notation N h = (S h + E h + I h + 
R represents the total human population contained in all 
the five compartments. We consider the solution of the 
system of equations (1) when the term 8 h \ h vanishes. In 
case if the death rate of humans due to malaria disease is 
considered to be free, i.e., 8 h — 0 then we 
obtain(d N h / dt) = (if/— pN h ). The solution of this 
differential equation is found to be N h (t) = (y//p) + 
[N h0 ~ (l^/p)] e ~ flt showing that N h (t) —> y//fi as t -» oo. 
The termAi h0 denotes the initial total human population. It 
can be interpreted that the total human population grows 
and asymptotically converges to a positive quantity given 
by ( y//fi ) under the condition that humans do not die due 
to malaria infection. Thus y//p is an upper bound of the 
total human population^ (t ) i.e.N h (co ) < i///fi. Whenever 
the initial human population starts off low below (y//p) 
then it grows over time and finally reaches the upper 
asymptotic value (y//p). Similarly, whenever the initial 
human population starts off high above (y//p) then it 
decays over time and finally reaches the lower asymptotic 
value (y//p) [3], 

Similarly on summing up all the individual equations from 
the system (1), it is straight forward to get dN m /dt — p — 
a) N m . Here the notation N m = (S m + E m + I m ) represents 
the total mosquito population contained in all the two 
compartments. The solution of this differential equation is 
found to be N m (t ) = (p/ca) + [N m0 — (p/w)] e~ Mt 
showing that N m (t) -> (p/o>)as t -> oo. The term 
Ai m0 denotes the initial total mosquito population. It can be 
interpreted that the total mosquito population grows and 
asymptotically converges to a positive quantity given by 
(p/w).Thus (p/ o>) is an upper bound of the total mosquito 
populationAi )n (t ) i.e. N m (oo) < (p/o>) . Whenever the 
initial mosquito population starts off low below (p/aa) 
then it grows over time and finally reaches the upper 
asymptotic value(p/w). Similarly, whenever the initial 


mosquito population starts off high above (p/oa) then it 
decays over time and finally reaches the lower asymptotic 
value (p/ca) [26], 

Hence all feasible solutions set of the human population 
and mosquito population of the model (1) enters the 
region. 

0 = i(S h , E h , 4 , R, S m , E m , I m ) £ Rl ; (S h , S m ) > 

0,Eh, /h, R,Em, Im>0,Nh<\f/ju,Nm<pa>. 

Therefore, the region /2 is positively invariant (i.e 
solutions remain positive for all times t) and in the model 
(1) is biologically meaningful and mathematically well- 
posed in the domain /2. 

3. Existence of Disease free equilibrium point E„ 

Disease free equilibrium points are steady state solutions 
when there is no malaria in the human population and 
there is no plasmodium parasite in the mosquito 
population. That is, absence of malaria causing infections 
occurs in both populations at the disease free equilibrium 
point. The disease free equilibrium point is denoted by 
E 0 = (Sh ,E h *,Ih ,R*, ,E m * Im). The equilibrium point is 
obtained on setting the right-hand side of the non-linear 
system (1) to zero. Thus, at the equilibrium point the 
quantities satisfy the condition E h * = 4 = R* = E m * — 
Im — 0, — (if//p) and Sm — (p/w).Also, Overhead star 

represents the values of the functions at the disease free 
equilibrium point. The disease free equilibrium point 
represents E^the disease free situation in which there is 
no malaria infection either in the society or in the 
environment. Therefore, the diseases free equilibrium 
point is given by 

E 0 = (y//p, 0,0,0, p/a>, 0,0) (2) 

4. Basic Reproduction Number R 0 

Reproduction number, denoted by R 0 , is the threshold or 
a level for many epidemiological models. It determines 
whether a disease can attack the population or not. The 
threshold quantity R 0 indicates the number of new 
infected individuals is produced by one infected 
individual. When R 0 < 0 each infected individual 
propagates the infection and produces on average less 
than one new infected individual so that the disease is 
expected to die out completely over time. On the other 
hand if R 0 > l,each individual produces more than one 
new infected individual so we would expect the disease to 
spread more and grow in the population. This means that 
the value of threshold quantity R 0 in order to eradicate 
the disease must be reduced by less than one. 

The following steps are followed to compute the basic 
reproduction number R 0 The basic reproduction number 
cannot be determined from the structure of the 
mathematical model alone, but depends on the definition 
of infected and uninfected compartments. Assuming that 
there are n compartments of which the first m 
compartments to infected individuals. That is the 
parameters may be vary compartment to compartment, 
but are the identical for all individuals within a given 
compartment. Let 

Xi = (x 1 , x 2 , x n ), X t > 0 for all, i = 1,2,. m 
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Be the vector of human and mosquito individuals in each 
compartment. Let us sort the compartments so that first 
m compartments infected individuals. 

Let F;(x) be the rate of appearance of new infections in 
compartment i. 

Vi(x ) = V~i(x ) — R + j(x) Where I/ + j(x) is rate of transfer 
of individuals into compartment i by all other means and 
V _ i(x) is the rate of transfer of individual out of the i th 
compartment. 

It is assumed that each function is continuously 
differentiable at least twice in each variable. The disease 
transmission model consists of non-negative initial 
conditions together with the following system of 
equations: 

V L (x), i = 1,2,3. n 


g = fi(x) = F l (x) 


dx ■ 

Where is the rate of change of x. The next is the 
computation of the square matrices F and V of 
order (mxm), where m is the number of infected classes, 


defined by F = 


r dFi f(x 0) ] and V = 


-dVi(x) 


(x 0 )] with 


dxj KKJ J J " L dxi 

1 < i,j < m, such that F is non- negative, V is non-singular 
matrix and x 0 is the disease-free equilibrium point (DFE). 
Since F is non-negative and V is non-singular, then V _1 is 
non-negative and also F V -1 is non-negative. Hence the of 
FV -1 is called the next generation matrix for the model. 
Finally the basic reproduction number R 0 is given by 
Ro = Y(FV- 1 ) 


Where y(A') denotes the spectral radius of matrix A and 
and the spectral radius is the biggest non-negative 
eigenvalue of the next generation matrix. Rewriting model 
system (1) starting with the infected compartments for 
both populations; S h , E h , l h , R, S m , E m , I m and then 
following by uninfected classes; S h ,R,S m also from the 
two populations, then the model system becomes 

Ut = — Y h - ^ + ^ E n 

dl h /dt = p h E h - (t + p + S h )l h 
dF / 

m Ut = Nh - (Pm + u)E m 
d\ m /dt S m co\ m (3) 

dS h /dt = y/+ cpR -—- pS h 


N h 


dR/dt = xl h — (<p + p)R 

, r / , _ 
uj m /dt — p ~ 


Since a h 


N h 


( s hm 0 I h'i 


a)S„ 


and a™ = 




malaria 


V N h J \ N h 

model. From the system of equation (3) F, and V t are 
defined as 




F(x) = 



V(x) = 


iPh + F)E h 

(t + \i+8 h )h -PhEh 
(Pm + C0)E m 
Pm 


The partial derivatives of (3) with respect to ( I h , / m ) and the jacobian matrix of F ; at the disease-free equilibrium point 
(2) is:- 

0 0 0 6 mh 0 I m 

0 0 0 0 
0/im0RP 


F = 


0 


a)Xp 


0 


0 


0 0 0 0 
Similarly, the partial derivatives of (3) with respect to ( E h , I h , E m , I m ) and the jacobian matrix Bps: 


v = 


(Pn + F) 0 

~Ph 

0 0 
0 0 


0 

(r + p + S h ) 
0 

~Pm 


0 


O) 


0 


v- 1 = 


(Ph + F) 


0 


(t + p + 8 h )(Ph + p) 
0 0 


0 


0 


(r + p + 8 h ) 


(Ph + F) 
Pm 

a>(Ph + a>) 


1 

O) 


FV- 1 = 


0 

Ph^h m 0PP 


MXp(r + \i + 8 h )(Ph +F) 
0 


0 

QhmQy P 
(r + p + 8 h ) 
Pn 


Pm Qmh® 

U)([] h + w) 


0 

0 


a> (Ph + w ) 


1 

0) 
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From FT -1 ,we can determine the eigenvalues of the basic reproduction number R 0 by taking the spectral radius 
(dominant eigenvalue) of the matrix FV _1 . Thus it is calculated by \A — AA\ = 0. We determine the expression for R 0 using 


the next generation matrix approach [18] as R 0 


P S m h9 hm 0 2 n 


I co 2 ^(r+A+n h + S h ) 

equilibrium point E 0 given by (4) is locally asymptotically stable if R 0 < 1 and unstable if R 0 > 1. 


Further, it can be verified that the disease free 


Local Stability of the Disease-Free Equilibrium 

The local stability of the disease-free equilibrium can be analyzed using the Jacobian matrix of the malaria model (1) at the 
disease free equilibrium point. Using [40], the following theorem contains 

Theorem:- The disease free equilibrium point for system (1) is locally asymptotically stable if R 0 < 1 and unstable if 
R 0 > 1 


Proof: The Jacobian matrix (/) of the malaria model (1) with 

S h = (E h + I h + R 8iS m = (E m + / m ) at the disease-free equilibrium point is given by: 

Ut = — T h - Wk + Wk 

dl h /dt = p h E h - (t + p + S h )\ h 
dR/dt = xl h — (cp + p)R 


dE m j _ Ih-bn 
'dt ~ Y h 

d\ m /dt S m a)l m 


(Pm + 9))E n 


(4) 


From the equation (4) the jacobian matrix 

~(Ph + AO 


(/) = 


Pn 


T 

a>xp 


0 


0 0 
- (r + p + 8 h ) 

~(<p + (o) 

0 

0 0 


'( Pm + w ) 
Prr 


Let a = + /r) , b — 8 mh 0, c — j3 h ,d = (r + p + 8 h ) ,e — t,/ = (<p + w ),g = 

(S m and j — a) 


g hmfiPP 

(j)1p 


9mhV 
0 

0 

0 

— 0 ) 

h= (Pm + 0>), i = 


Thus 


(/) = 


0 
L 0 


e 

9 


0 

0 


—h 


The eigenvalues of jacobian matrix are:- 

— (n + T) 


(/) = 


0 


e 

9 


17 - A/| = 0 
0 

d 0 


/ 


0 


0 


-(/i + A) 


b 1 


0 

(j + X) 


The third column has diagonal entry, therefore one of the eigenvalues of the jacobian matrix is —(cp + a>) or/. By using 
Routh-Hurwitz, stability criterion. 


Bn = 


P9 m h9hm ® P 
«V(t +^+V-h+ s h) 


the proof end. 


Using the Routh-Hurwitz criterion is a method for determining whether a linear system is stable or be examining the 
locations of the characteristic equation of the system. In fact, the method determines only if there are roots that put 
outside of the left half plane; it does not actually compute the roots Routh-Hurwitz criteria [12]. 


We can determine whether this system is stable or not, checking the following conditions:- Two necessary but not 
sufficient conditions that all the roots have negative real parts are 

> All the polynomial coefficients must have the same sign. 

> All the polynomial coefficients must be nonzero. 
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The necessary condition that all roots have negative real 
parts is that all the elements of the first column of the 
array have the same sign. The number of changes of sign 
equals the number of roots with positive real parts. All the 
elements of a particular row are zero. In this case, some of 
the roots of the polynomial are located symmetrically 
about the origin of the A plane, e.g.,a pair of purely 
imaginary roots. The zero row will always occur in a row 
associated with an odd power of A. The row just above the 
zero row holds the coefficients of the auxiliary polynomial. 
The roots of the auxiliary polynomial are the 
symmetrically placed roots. Be careful to remember that 
the coefficients in the array skip powers of A from one 
coefficient to the [2], 

5. Existence of Backward Bifurcation 

We intend to determine the stability of the endemic 
equilibrium and to carry out the possibility of the 
existence of backward bifurcation due to existence of 
multiple equilibrium and reinfection. As a disease attacks 
it reduces the number of susceptible individuals in the 
population, which tends to reduce its reproductive rate. 
For a backward bifurcation to occur, this means that when 
R 0 < 0 the endemic equilibrium point can exist as well as 
when R 0 > 1. 

we would expect the disease to be able to attack at R 0 = 
1 in the case of a backward bifurcation with the properties 
of unstable equilibrium bifurcating from the disease-free 
equilibrium when R 0 < 1,giving rise to multiple stable 
states. But not in the case of a forward bifurcation, in 
which in the absence of a low-level unstable equilibrium 
when R 0 < 1 and a stable equilibrium bifurcating from the 
disease-free equilibrium when R 0 > 1 arise naturally 
when the disease does not attack when R 0 = 1. A simple 
criterion for a backward bifurcation, then, is one in which 
the disease can attack when R 0 = 1. This implies that the 
disease-free equilibrium may not be globally 
asymptotically stable even if R 0 < 1. 

6. Numerical simulation 

In this section we consider the simulation study of the 
system of differential equations given in [1], As stated 
earlier these equations describe the dynamics of human 
and mosquito populations of the malaria model that 
includes intervention strategies. The simulation study is 
performed using ode45 solver of MATLAB software. The 
Runge - Kutta fourth-order method based on a variable 
step-size is used for the purpose. The parametric values 
have been collected from the literature and used here. 
Those were not available were not obtained from 
literatures published by researchers in malaria endemic 
countries which have similar environmental conditions, 
we present the numerical analysis of the model(l). The 
initial conditions used were S h (0) = 47186 ,E h (0) = 
16987,4(0) = 47473,£(0) = 47470, S m (0) = 

17500,£4(0) = 8750,/ M (0) = 26,250. We simulate the 
basic malaria model in the absence of any intervention and 
the malaria model without intervention strategies, and 
find out the effects of varying each intervention parameter 
[60],[59], 

6.1. Estimation of Parameters 

We estimate that it will take 3 times a day for 7 days to 
recovery from malaria infection through Chemotherapy 


and the incubation period of malaria in humans was 
considered [9],[39]and [59], 


Table 3: Estimated Parameter Value of the Malaria 
Modal without Intervention Strategies and with 
Intervention 


Symbol 

Values 

sources 

xp 

0.000027 

Calculated 

@hm 

0.0055 

[60] 

P 

0.071 

Niger,2008 [24] 

®mh 

0.42 

Estimated 

Ph 

1/14 

FMOH.2004 [9] 

<P 

0.40 

chitnis 2005 [58] 

An 

1/11 

chitnis 2008 [24] 

A 

0.01 

Miranda. 2009 [30] 

T 

1/7 

FMOH.2004 [9] and Tamwiine et al 2004 [18] 

€ 

0.00722 

Gumel,2008 [12] 

P 

0.0000548 

Calculated 

g 

0.78 

Mwamtobe,2010 [34] 

LJ 

1/25 

Blayneh,2009 [4] 

v 

0.457 

[59] and [4G] 

6 

0.00000071 

World malaria report 2014 [45] 

7 

0.11 

Mwamtobe,2010 [34] 

‘P 

1/121.75 

Estimated 


6.2. The System of Human Population State Variables 
of the Basic Model without Intervention 

The simulation of basic model has been conducted to find 
out the dynamics of the disease in the population when 
there is no intervention to reduce or eradicate the disease. 
In the absence of interventions strategies, the susceptible 
populations in Figure 2 as shown below red colored that 
the change in state variables of malaria model shows the 
dynamics with time of susceptible humans with 
R 0 = 1.3874. 


Plot of human population against time 



Figure 2:- A phase portrait illustrating the changes in 
the four state variables of the malaria model showing 
the system with time, of susceptible humans, exposed 
humans, infected humans and shows the system of 
recovered humans with R 0 — 1. 3874. 
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This susceptible population increases then remained 
constant within the three years. The Figure 2 were blue 
colored showing the system with time of exposed humans 
initially was high but progress it reduces as some of 
population enters infections class and other recover. The 
green colored also were indicated in the Figure 2 showing 
the dynamics with time of infectious humans the infection 
population reduced as some recovers and other die within 
this three years. Finally, the yellowed color in Figure 2 the 
dynamics with time of recovered humans the recovered 
population initially was high as time goes by the recovered 
people reduced or to goes back to susceptible class. 



Time(years) 


Figure3:-Illustrates the changes in the three state 
variables of the malaria model showing the dynamics 
with time, of susceptible mosquitoes, exposed 
mosquitoes and infectious mosquitoes with 
R 0 = 1.3874. 

In Figure 3 shows all the three curves are decreasing as 
time increases which are positive for the current 
interventions in the mosquito population, but there is still 
more work to be done in the human population. Therefore, 
we will consider the effects of varying the main 
parameters responsible control malaria after considering 
malaria prevalence rate in the population now. 

Prevalence the Basic malaria model without 
Intervention Strategies 

Prevalence is defined as the ratio of which the number of 
cases of a disease in a population and with the number of 
individuals in a population at a given time. 



Figure 4:- Represents changes of prevalence with time 
with R 0 = 1. 3874. 

The prevalence graph shows that the prevalence rate as of 
now is high which confirms the Figure 2 as shown above 


that there is more work to be done if we want to achieve 
malaria free society, because the prevalence rate reduces 
asymptotically to zero as possible to increasing the 
intervention strategies time to time. 

6.3. Simulation of Biting Rate of Mosquitoes on the 
Basic malaria model without Intervention 
Strategies 

We now consider the effects of varying the main 
parameters responsible for controlling the spread of 
malaria disease. The values of the biting rate of 
mosquitoes, transmission rate of infection from an 
infectious mosquito to a susceptible human, rate of loss of 
immunity for humans and the mosquito population were 
reduced by constant fraction 1/8 , while the values of the 
other parameters are maintained. This important to 
showing the relationship of the susceptible and infected 
human populations was considered. 



Susceptible individuals 


Figure 5:-Illustrates the system of infected human 
population against susceptible human population with 
i? 0 =0.1733 

The Figure 5 with R 0 =0.1733 shows that at an initial 
stage, the susceptible human population was free from the 
disease. The infection in the population is shown that it 
increases with time when there is no intervention being 
practiced. As most of the susceptible individuals were 
getting infected with time, the infected individuals 
increases. This is evidenced further when the parameter of 
mosquito biting rate varied. 

7. Conclusion and Recommendation 
7.1. Conclusion 

Analysis of the model showed that there exists a domain 
where the model is epidenriologically and mathematically 
well-posed. The important parameter in our model, the 
basic reproduction number R 0 as an improved control 
intervention measure was computed. The model was then 
qualitatively analyzed for the existence and stability of 
their associated equilibria. It was proved that under the 
condition that R 0 < 1 the disease-free equilibrium E 0 is 
locally asymptotically stable, and when R 0 < 1 the 
endemic equilibrium E 1 appeared. The model exhibits the 
phenomenon of backward bifurcation where a stable 
disease-free equilibrium co-exists with a stable endemic 
equilibrium for a certain range of associated reproduction 
number less than one. 
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The numerical analysis of the model suggested that the 
most effective strategies for controlling or eradicating the 
spread of malaria disease. The use of insecticide-treated 
bed nets and indoor residual spraying and prompt and 
effective diagnosis and treatment of infected individuals. 
This study agree [7] suggestion that the intervention using 
insecticide-treated bed nets represents an excellent 
example of implementing an infectious disease control 
programme, and [38] study, which showed that both 
regular and non-fixed spraying resulted in a significant 
reduction in the overall number of mosquitoes, as well as 
the number of malaria case in humans. 

7.2. Recommendation 

So, far no really effective vaccine has been developed 
against malaria, so we cannot protect ourselves against 
the disease. For many years there have been few effective 
treatments for malaria, but things are getting better. In a 
country like ours, drugs alone are the answer. Many 
peoples live a long way from medical centers and cannot 
reach them easily, and medicines can be expensive. 
Mathematical modeling of spread of malaria disease can 
provide understanding of underlying techniques or 
strategies for the disease spread, help to pinpoint key 
factors in the disease transmission process, suggest 
effective controlling and prevention measures, and 
provide estimate for the severity and potential scale of the 
endemic and epidemic. The following recommendations 
should be considered: 

> Individuals should be award about severity of malaria 
and increasing personal protection measures are 
highly recommended to prevent malaria [to control 
the reproduction number R 0 < 1 ]. 

> using methods of controlling malaria must involve 
controlling the Anopheles mosquitoes. Whenever 
possible avoid contact with mosquitoes. 

> Using mosquito repellents, having screens on doors, 
windows to prevent mosquitoes, wearing clothes that 
protect the skin against mosquito’s long sleeves and 
trousers are all effective measures that can protect 
you against malaria. 

> Well-made insecticide-treated bed nets and indoor 
residual spraying and prompt and effective diagnosis 
make a big difference and they are cheap treatment 
effective way to control malaria transmission. 

> Minimize any opportunities for the mosquito breed. 
They will lay eggs in any standing water in garden 
pond, old trees, flower pots, old drink cans,...etc. 

> Removing the mosquitoes breeding places by 
removing as much as standing water as possible. The 
simplest way to do this to make sure you store 
rubbish out of the rain and dispose of your rubbish 
properly. 

> Proper disposal of sewage again, managing human 
waste so that foul water is not left around will reduce 
the breeding places for the mosquitoes. 

> Biological control [where an organism that feeds on 
the larva is introduced to the water] and chemical 
control [pesticides or IRS] spraying on to the water 
where the mosquitoes breed will kill the eggs and 
larvae, this in turn reduces the number of mosquitoes 
and slower the infection rate. 

> Because of the complications of measuring malaria at 
different levels with different immunological status 
prevalent in different age and gender groups, and 


across different locations, some guidelines should be 
developed to give researchers and health 
professionals a more accurate foundation on which to 
select their indicators. 

> There should be administration of drugs to these 
malaria case encountered individuals irrespective of 
whether symptoms show up or not in order to 
decrease the natural birth rate of mosquito and 
increase the rate mosquito using Indoor-Residuals 
Sprays [IRS], 
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